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INTRODUCTION 

A study of the effects of earthquakes on the rotational motion of the earth has 
been conducted. The analytical developments providing the connection between 
the fault parameters and the corresponding changes in the moments and products 
of inertia are due to Rice and Chinnery (1972), the method involves the applica- 
tion of the reciprocal theorem of elasticity and Volterra's formula as well as the 
displacement and stress fields for the second degree static response of the earth 
model being used. 

Two earth models have been used in the investigation, the parametric model 
due to Dziewonski, Hales and Lapwood (1975) and the M 3 model of Landisman, 
Sato and Nafe (1965) as given by Israel, Ben-Menahem and Singh (1973). 

In order to obtain the displacement and stress fields it is necessary to inte- 
grate numerically a system of differential equations representing the state of 
equilibrium of an elastic body. The numerical integration problem presents cer- 
tain aspects which require consideration; some of these aspects are the following: 
the conditions in a neighborhood close to the origin, at which the differential , 
equations become singular; and the appropriate boundary conditions between the 
solid and the liquid parts of the body. 

The numerical results of the investigation yield the magnitude and direction 

of the pole shift as well as the change in the length of the day. In addition, the 

* 

changes in the second degree coefficients of the geopotential have been .computed. 
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In order to generate numerical results the source parameters corresponding 
to the Alaskan earthquake on March 28, 196,4 were chosen, as given by Israel 
et al. , 1973. 


1. RELATION BETWEEN' MOMENT OF INERTIA, ROTATION AND FAULTING 
Consider a body of mass M rotating with angular velocity co and let the axis 
L be defined by the direction of w. 


The moment of inertia with respect to the axis L is, defined by 


■/ 


I L = / d 0 dm 

M 


( 1 . 1 ) 


d 0 is the perpendicular distance from the element of mass dm to the axis L. 

* 

Assume a displacement field to be defined throughout the body due to causes other 
than rotation, i.e. , as a result of faulting. The perpendicular distance from dm 
to the axis L is now given by 

d = d 0 +Ad (1.2) 

where Ad denotes the component of displacement due to faulting perpendicular to 
L. The moment of inertia is now given by 


I 


L 



(1.3) 


substituting Equation (1, 2) into Equation (1. 3) yields, 


1 L = / d^dm+2 I (d 0 Ad)dm + / Ad 2 dm 

•'m •'KA 


(1.4) 
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Assume Ad to be small enough so that the term containing Ad 2 can be neglected-, 
then 


AI l = 2 


(d 0 Ad)dm 


(1.5) 


where AI L stands for the increment in the moment of inertia due to faulting. 
Now, the rotational potential is given by 



( 1 . 6 ) 


and the corresponding force by 



co 2 d 0 


(1.7) 


Therefore, 


d o 


_L 1L 

CO 2 3d fl 


substituting Equation (1.8) into Equation (1.5) yields 


AI 


L 



dm 


( 1 . 8 ) 


(1.9) 


2. RECIPROCAL THEOREM AND VOLTERRA'S FORMULA 

The Reciprocal Theorem of elasticity is due to E, Betti, its proof can be 

* \ 

found in the book by Love (1944), more recently Smylie and Mansinha (1971) and 
* 

Rice and Chinnery (1972) have extended it to the case of self-gravitating elastic 
systems in states of large initial hydrostatic stress. The presentation which fol- 
lows is due to Rice and Chinnery (1972). 
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Consider a body of mass M and a surface £ within the body, define the 
following: 

n + , IT : unit vectors normal to an element of surface d2, opposite in direc- 
tion, with n _ pointing from d£~ to d£ + , 

IT ; a displacement field defined throughout M with a discontinuity, 

Au =tt + -u~, due to slip, defined on 2. 

[o-] : stress field defined throughout M, due to* tT. 

f* : a body force field, defined throughout M, independent of u. 

u* :• a displacement field defined throughout IVI, due to the application of 

f*. 

[a*l : stress field defined throughout M, due to u*. 

Assume the following conditions to holdj 

i. [a], and [o*] vanish at the outer surface of the body.. 

ii. [o*]. and u* are continuous across. 2, since they are produced by f*, with 
no slip on 2. 

iii. (n, • [a] ) + + (n * [u] )“ = 0, this last condition implies, equilibrium after 
slip. 
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The Reciprocal Theorem states that the work done by the body force f* and 
the stress [a*] acting through the displacements tT is equal to the work done by 
by the stresses [a] acting through the displacements tT*. Therefore, 



; ( 2 , 1 ) 


By virtue of conditions (ii, ) and (iii.) Equation (2.1) can be written as: 


/ ( 7 * * u)dm = / (n • [o*] * Au)d£ 
•'2 


( 2 . 2 ) 


with n = "n , now let Au" = (Auys, s' being a unit vector in the direction of slip. 
Furthermore define 


r* - ft • [o*] • t 


(2.3). 


Equation (2. 2) can then be written as 


I (f* * u)dm = j r*(Au)d2 
•'m •'s 


(2.4) 


Equation (2.4) is referred to as Volterra’s formula. 


3. RELATION BETWEEN MOMENT OF INERTIA AND VOLTERRA'S FORMULA 
Consider the integrand in the left side of Equation (2.4), let f* denote a body 
force due to rotation as given by Equation (1. 7) and let u stand for the displace- 
ment field due to faulting. Equations (1. 9) and (2. 4) then yield 


•5 
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where 


AI l = ~ / r*(Au)dS 


r£ = 11 * l>*] L *'s 


(3.1) 


(3.2) 


[a*] L being the stresses due to the rotation about the axis L. The Stress tensor 
[a*j L contains a factor oo 2 which will cancel the same factor appearing in the 
denominator on the right side of Equation (3* 1), this will be shown below. 


4. MOMENTS AND PRODUCTS OF INERTIA 
Consider the inertia dyadic. 



r * 

W 1 

I ^ 

~W } 

-i, z tr 

[I] * 

/s>/\ 

-v 1 

IyyJ'J 



T 

~^xz^ 1 

“^2^ J 



Note that. 


f- fu; m = L 


k - m • k =• L, 


Define the following unit vectors: 


rn : =' Wrf 3$^ 


£ =' 0 


n = (d'A/2)(t-j)- 


(4, i) 


(4.2)' 


(4.3)- 


6 - 
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Then, 


q = (iA/2)(1c+1) 


-m • [I] -m-r • [I] -r“ = 2L 


P*[I] -p-t*[I] ;t= 21y 2 
o - [I] *o-q • [I] • q 5 s -21 


Equations (3. 1), (4. 2) and (4. 4) then yield 


AI XX - — I ^(Au>dS 

“ •'S 


A, yy ‘ -i / 


AI zz “ ~ / r*(Au)d2 


f t^(A u)dS-— f r* 

•'y ^0“ y y 


AI xy =" J ^(Au)dS- 
co ^2 


(Au)dS 


AL= 4r / ^(Au)d2 rf(Au)dZ 

J co" -^2 F co ^2 




L 1 * 


As indicated by Equation (3.2) the quantities r*, r*, . , r* denote the shear 
stresses induced on 2 by rotation about the axis indicated by the subscript. 


If the seismic event is assumed to be represented by a point source then 


Equations (4.5) can be written as follows. 



10 


where 


AI XX = 2K r* 

AI yy - 2K f* 

AI ZZ - 2kr* 

AI X y = K(r* -r*) 
AI yz = K(r*-r*) 
AI XZ = K(t* -t*) 

K = (1/w 2 )(AuX2) 


(4.6) 


'(4.7) 


5. STBESSES 

Equations (4. 5) and (3. 2) require knowledge of the stresses produced on 2 
by rotation about a certain axis, this constitutes a part of a more general prob- 
lem which, is. briefly outlined below.. 



U = (u x ,u y ,u z ) (5.3) 

The body force potential V is given by 

V = V 0 + * (5.4) 


8 . 
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where denotes the perturbation from the- "undisturbed state V 0 . V must obey 
Poisson’s equation: 

V 2 V = -4vrGp (5.5) 

Let the disturbing potential be given 

T n = (r/a) n S n «?,tf) (5.6) 

4 

where S n denotes a surface spherical harmonic and a the radius of the earth. 


Equations (5. 1) and (5. 5) can be expressed in spherical coordinates and with 
the appropriate set of boundary conditions they will admit the following eigen 
vector solutions: 


u r = y,CDS n 


dS n 

u » = v ’ (r) Y 


y 3 w 3s„ 

^ - sin0 d(j> 


u* = 


(5.7) 


The stress-strain relations yield: 

*rr = y 2 W S n 
3S„ 


a x6 = y 4 w 


de 


°r<p ~ 


y 4 fr) 9Sn 

sin 9 d<f> 


(5.8) 


f2 n(n 

a eo 55 [7 ( x + ^)yi + ^y\ ~ 


i+l) , 1 c ^ /2p \ 3 2 S n 

r V3 J " l r y3 J O0 2 
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where. 




2 . n(n 

— (X Ay.j- 




^ * 4 '* 

+ cotey ys j-* y y^ 

\ r V * 30 \rsiri 2 0 7 30 2 

_ / 2 m cotO \ 9S„ / J* 

\ r sin 0 y */ d<p \rsin 2 '0 7 


d<p 2 


* (r^TTT ^ 


9030 


y 3 


= X j*, +^- 


n(n + I ) * 

yj — — y 3 


+ 2juy, 


(5.8) 
(cont. ) 


(5.9) 


- (■ y 3 + y i \ 

y 4 (y 3 - T+T ) 

The time-dependent solution e lfnt has been omitted from Equations (5,7) and 
(5.8). The dots denote (d/dr). The stresses given by Equation (5.8) are those 
produced by the action of the disturbing potential T n and do not include the initial 
hydrostatic stress. 


6. BADIAL FUNCTIONS 

The radial functions y t * y 2 , y 3 and y 4 appearing in Equations (5. 7) and (5.8) 
have to be obtained by numerical integration when working with models resem- 
bling the real Earth. In addition define 


0 = y 5 (r)S n (0, d) 

(6.1) 

y<, = y 5 -4*<<yi 

(6.2) 


10 



Equations (5.1) and (5.5) can then be expressed as a system of six ordinary 
differential equations (Alterman et al. , 1959) 


13 


t = [M]y 

(6.3) 

y = (y 1 ,y 2> y 3 ,y 4 ,y 5 ,y 6 ) T 

(6.4) 


The static or steady-state solution is obtained by letting the frequency f n 
be equal to zero. The elements of the matrix M are given by Smylie and 
Mansinha (1971) and Israel et al. (1973). They are given below for the' sake of 
completeness. The symbols p Q , p and X denote the density and the elastic param- 
eters, g 0 stands for the gravitational acceleration and n is the degree of the de- 
formation, in this case n = 2. 

Let 


c = l/(\ + 2p) 

l - 

M u = -2Ac/r 
M 12 = c 

Mj 3 = n(n + 1 )c/r 

M 71 = 4p(3\ + 2fi)c/r 2 -4p 0 g 0 /r 

(6.5) 

M 22 = -4pc/r 

M 23 = -n(n + 1 )[2/x(3A + 2p)c/r 2 - P 0 g 0 /r] 

M 24 = n(n+ l)/r 

M 26 ~ p o 

M 3 , = -1/r 

M 33 = 1/r 

11 
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M 34 = 1/a 

M 4 ‘ = -2ju(3X + 2n)clx 2 + A 0 g o /r . 

M 42 = -Xc/r 

M 43 = 4n(n+ l)p(X + p)c/r 2 - 2p/r 2 
M 44 = -3 /t 

m - / (6 * 5) 

^45 Pj x (cont.) 

M 5i = 4nGp 0 

M 56 = 1 

M 63 = -47rGn(n + 1 )p 0 /r- 
M 6S = n(n + 1 )/f 2 
M 66 - -2/r 

all the other elements are equal to Zero. 

Equations (6. 3) and- (6. 5), are applicable to the solid, regions of the, earth. 
The equations for the liquid core have been developed, by Smylie, and Mansinha 
(1971) and Israel et al. (1973), they are given below: 

yi = y 5 /g 0 

y 2 = ° 

y 3 = (4y s + ry 6 )/n(n i + i)g 0 ^ 6 * 6 ^ 

y 4 -o 

y 5 = (47rGp 0 y 5 )/g o + y 6 

y 6 = [n(n + 1 )/ r 2 - 16irGpJg q T]y 5 -(47iGp Q /g Q + 2/r>)y 6 


12 
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The following assumptions are implicit in the set of Equation (6.6): 

i. the liquid core is in a state of hydrostatic equilibrium before and after 
the deformation, this means that the tangential stresses are equal to 
zero, 

ii, the dilatation and the normal stress are equal to zero. Zero dilatation 

is consistent with constant core volume and makes the Adams-Williamson 
condition unnecessary. 

Equations (6.3) and (6.6) have to satisfy certain conditions at the boundaries be- 
tween the solid and the liquid parts of the earth, there has been some debate in 

> 

the literature concerning this issue. The conditions given below are those de- 
veloped by Israel et al. (1973) and Crossley and Gubbins (1975), at least with 
respect to the condition for the variable y 6 . Let {yj} denote the jump (discon- 
tinuity) in the variable yj at the boundary between the liquid and the solid parts, 
then 

{ y i} = some constant 

{^2} = PoZo frl} 

|y 3 j. = some constant (6.7) 

W = 0 
{y 5 } = o 

{yj = -4 ttG Po {y,.} 


13 
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Note that Equations, (, 6. 7)fimply the continuity* oPy 4 ’ an dy 5 . ' 


Alterman et al. (1959) developed the remaining, boundary conditions to be 
imposed on the differential 1 equations. At. the- deformed surface- of the earth the 
stresses must vanish and the gravitational’ potential, and its: gradient must be 
continuous, i.e. , 


y, 2 (a) - 0 
y 4 (a) - 0 


y 6 (a) + 


(n + 1 ) 
a 



b. 


( 6 . 8 ). 


where a denotes the. radius of the earth and b is; a. constant witl*- a value depending 
on the nature of the disturbance as. shown by Takeuchi et at, (i9.62) r In particular^ 


for the earth-tide problem: 

b - (2n+-l)/a (6. 9) 

Implicit in Equation (6. 9) is the assumption that. the, disturbing potential is givep 


by Equation (5.6). 


7. THE ROTATIONAL POTENTIAL 

The rotational potential as giyep by Equation, (l t 6), gap- be pjqmessed in the. 
form of Equation (5.,6_) (Sanchez, 19,74): 

T 2 -- (7.1) 

where 


2 

S,(0, <j>) = ^ P™ (cos 6 )(q™ cos m$ +<- u^ 1 sin jnb) (7.2) 

m~= o 
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q° = (a 2 /6)(co£ + co 2 -2 co 2 ) 
qj = ~(a 2 /3)(co x co z ) 

u 2 = ~C a2 /3)(<^yW z ) (7.3) 

q 2 = (a 2 /12)(co 2 -coj) 

u 2 = -(a 2 /6)(co x co y ) 

Equation (7.1) does not include a term equal to (1/3) (co 2 + to 2 + co 2 ) which pro- 
duces expansion or contraction of the body as a whole, such a deformation will 
affect the moments of inertia. 

The stress tensor which appears in Equations (3. 2) and (4. 5) requires the 
evaluation of the stresses given by Equation (5. 8), the stresses are produced by 
rotation about a certain axis, the components of rotation co x , w y and co z appear- 
ing in Equations (7.3) will take values corresponding to the axis: 

(co x , co y , w z ) x = (co, 0, 0) 

(co x , to y , GJ z ) y = (0, co, 0) 

(co x , co y , co z ) z = (0, 0, co) 

(co x , co y , co z ) r = (lA/2)(w, co, 0) ^ 

(co x , co y , co z ) m< = ( 1 A/2) (co, -w, 0) 

(co x , C 0 y, co z ) g “ (lA/2)(0, co, co), 

(co x , co y , co z ) p = (lA/2)(0, -co, co) 

(co x , co y , co z ) 0 = (1 A/2)(-co, 0,,co) _ 

(co x , co y , C0 7 ) q - ( 1 />/ 2)(co, 0, co) 


15 
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8. SOURCE PARAMETERS 

» t 

The shear stresses given by Equation (3. 2)' contain the source parameters 
n and s' as well as the stress tensor due to rotation. The stresses given by 
Equations (5.8) are computed in a spherical system of coordinates, the source 
parameters n and 's are usually given in an epicentral coordinate system, the 
necessary transformations are given below: 
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5: dip angle 

X: slip angle 

a: azimuth angle 

9. OUTLINE OP PROCEDURE 

The numerical results which are given in the last section were obtained by 
assuming the seismic event to be represented by a point source. The following 
is an outline of the procedure. 

A. Given h (source depth) integrate the -system iof (differential ^equations 
given by Equations - s(6,.6) -to obtain '(y^ „ y 2 ., y 3J ly-j) ;at the cor- 
responding value of r. 

B. Choose an axis of rotation, the corresponding set.'(;oo x ., co y ., co z ) will be 
given by Equations (7.4)., ,co yean be ;set -equal -to -1, 

C. Using the results of step B and Equation >(.7. 3) compute Q 2 * u l> 
u 2 for each axis of rotation. 

D. Given 6 and <t> (col'atitude and longitude of the source) tcompute S 2 and 
its derivatives for each rotation axis by means of Equation >(#.',2) and the 
results of step*C. 

E. Using the results of step's- A and D and Equations. (5. 8); compute a rr , a^,, 

°08’ a W' °8<!> :for each station axis. 


: 18 



20 


F. Given 5, a (dip, slip and azimuth of the fault) compute s r , s 0 , s^, 
n r , n e , as given by Equations (8.4). 

G. Using the results of steps E and F and Equation (8. 1) compute t* for 
each rotation axis. 

H. Given Au, E (slip magnitude and fault area) and using the results of step 
G compute the changes in the moments and products of inertia by means 
of Equations (4.6). 

I. Use the results of step H to compute the magnitude and direction of the 
pole shift as well as the change in the length of the day. It is possible 
also to compute the changes in the second degree coefficients of the 
geopotential. 

10. NUMERICAL PROCEDURE 

The numerical integration of the system of differential equations given by 
Equation (6.3) - (6.6) requires the adoption of an earth model giving the radial 
distribution of density p and the elastic parameters ju and A. The two earth models 
used in this investigation are the following: 

i. the parametric earth model due to Dziewonski, Hales and Lapwood 
(1975) in which radial variations of density and elastic parameters are 
represented by piecewise continuous analytical functions of the radius. 
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In particular the model representing: a- continental structure was used 
and it is referred to’ as tlie P model, below. 

ii. the model M 3 of Landisman, Sato 1 and- Nafe 1 (.1-965) as given by Israel, 
Ben-Menahem and Singh (1973)’. The density and elastic parameters- for 
this model are given in table form and a cubic spline interpolation was 
used to obtain their values as- functions of the radius. 

The numerical integration of model P is performed 1 as : follows: 

a. a homogeneous earth model is used' to obtain 1 nominal values- for the. 

radial’ functions (y 2 , y 2 , y 3 , y 4 , y 5 , y 6 ) at a point within the inner - solid’ 
core,, i. e. , r 0 = 6..371 x 10 6 cm. The. system of differential'- Equations. 
(6. 3) is numerically integrated' outwards 1 to the boundary between 1 the- 
solid inner core and the liquid' core., Sinee y 2 vanishes! in? the. liquid 
core. Equations (6.7). yield:. 

{y 2 } - -y 2 - 

{yj = (:l/p 0 .g b > {y 2 }, 

(I'Owl)- 

K}. = M 
W - °' 

The third Equation (6..6); yields, the value o£y 3; w.ithin.the- liquid: core. 
There- remains- the condition matching the* values of’y 4 , at the boundary,, 
i.e.,_ 
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b. the homogeneous earth model yields nominal values for the variables 
(y t * y 3 1 y 5 ) at the surface of the earth where y 2 and y 4 are required to 
vanish. The value of y 6 is obtained from Equations (6.8) and (6.9), 
Equations (6.3) are then integrated inwards to the liquid core boundary. 


e. the two sets of solutions meeting at the liquid core -mantle boundary 
must satisfy the conditions given by Equations (6.7), i.e., 


=0 

(y 4 ) = 0 

LC ' . (10.3) 

LC^M = ° 

=0 


d. a general purpose adaptive iterator for nonlinear problems (Campbell 
et al. ) is used in order to satisfy the five conditions expressed by Equa- 
tions (10.2) and (10.3), which play the role of dependent variables. The 
independent variables subject to variation are the values of (y ^ , y 2 , y 3 , 
y 4 , y 5 , y 6 ) at r 0 = 6. 371 x 10° cm and the values of (y, , y 3 , y 5 ) at the 
surface of the earth. The numerical integration package consists of an 
Adams-Moulton, Runge-Kutta fourth order. 


21 



23 


The numerical integration of model M 3 follows a similar procedure but there 
are some differences which are outlined below: 

a. the absence of a solid inner core in this model requires that nominal 
values of (y 5 , y 6 ) obtained from the homogeneous solution at r Q = 6.371 

x 10 6 cm be used in the integration of Equations (6.6) throughout the liquid 
core. 

b. the value of the gravitational force, g 0 , is obtained from the integration 

of an additional differential equation (Pekeris and Jaroseh, 1958) i. e. , 

dg 0 2 

-2- + - g 0 = 4 ttGp 0 (10.4) 

dr r 

Nominal initial values of g 0 for the integration of Equation (10. 4) are ob- 
tained from the homogeneous solution at r Q = 6. 371 x 10 6 era and at the 
surface of the earth. Model P does not require this procedure since the 
density is given by polynomials, 

i 

e, the adaptive iterator uses the values of (y 5 , , g Q ) at r Q and the values 

of (y ,, y 3 , y 5 , g Q ) at the surface of the earth as independent variables 
in order to satisfy the set of conditions expressed by Equations (10.3) and 
the condition on g 0 at the liquid core-mantle boundary, i. e. , 

L C fe}„ - 0 <10.5) 

Once the adaptive iterator has converged to a set of initial conditions the so- 
lution trajectories for the radial functions have been obtained. Figures 10. 1-10. 12 
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below give the results for the P and M 3 models^- The units for the radial func- 
tions are omitted in the figures, they are given by Equations (5.7), (5.8), (6.1) 
and (6. 2). 

11. NUMERICAL RESULTS AND CONCLUSIONS 

The magnitude of the pole shift can be obtained from the following equation 
(Ben-Menahem and Israel, 1970): 

(ICOt; . 

MS = — - (AI* Z + AI^)* (11. 1) 

Af o • • 

where a is the radius of the earth, <o E the -mean angular velocity of the diurnal 
rotation, f Q is the angular frequency of the free Chandler wobble and A is the 
mean equatorial moment of inertia. The following' values were used in the 
computation 

a - 6.371 x 10 8 cm 
w E = 7.292 x 10“ 5 rad/see 
f Q = 1,69 x 10“ 7 rad/sec 
A = 8. 016 x 10 44 gm cm 2 


The direction of the shift is given by 


AS = arc tan 



.( 11 . 2 ) 


The change in the length of the day is calculated by means of the following equa- 
tion (Munk and MacDonald, 1975, page 98) 
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AL V 

ALOD =* (LOD) . (lli 3) 

C • f - • 

where LOD is the ‘mean length of the day and & is the mean polar moment of in- 
ertia. The following values were used in the computation: 


LOD = 86,460 sec 
G = 8.043 10 44 gmcin 2 


The changes in the. second degree coefficients of the geopotential are given by well 
* 

known relations which are reproduced Below: 

AC? = (Q/a)(Al x . x tAi yV '-2AI zz ) 

AC* = Q(Al k2 ) 

ASj = Q,(A^ 2 ) 

( 11 . 4 ) 

AC? = (Q/4)(AIyy - Ai xk ) 

AS? - (Q/2)(AI ky ) 

Q = l/a 2 M 

Where M is the mass of the Earth* 

M - S.975 x 10 27 gm » 

In order to generate numerical results a seismic event has td be chosem 
The Source parameters corresponding to the Alaskan earthquake oh March 28* 
1964 are used, as given by Israel et al* * 1973* the parameters giving the po- 
sition of the source are the followings 


36 



26 


6 = 28.9° 

<{> = 212.4° 
a = 225° 

where d denotes colatitude, 0 the longitude and a the fault azimuth. The depth of 
the source is allowed to take the values of 20, 60, 100 and 200km, and numerical 
results are given for each ease. The magnitude of the source is determined by 
the slip magnitude Au and the fault area 2, 

Au = 20 meters 
2 = 70,000 km i. 2 

The source mechanism is specified by the dip angle 5 and the slip angle A. Nu- 
merical results are generated for the three assumptions below, 

i. vertical strike-slip: A = 0°, 5 =90° 

ii. vertical dip-slip: A = 90°, 6 = 90° 

iii. dip-slip on 45° plane: A = 90°, 5 = 45° 

Tables 11. 1-11. 8 below give the results generated for the parametric earth 
model due to Dziewonski et al. denoted as model P, and the earth model M 3 of 
Landisman et al. 

The results indicate that the magnitude of the pole displacement depends on 
the assumptions concerning the mechanism of the source, the same can be said 
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With respect to the direction of the shift and the change in the length of the day. 

The type of earth model used in the computations has some bearing in the results 
as is to be expected but not nearly as much "as the type of source mechanism. The 
depth of the source introduces variations in the magnitude of the displacement and 
the change in the length of the day but the direction of shift is insensible to this 
parameter in the cases corresponding to a vertical strike-slip and a vertical dip- 
slip, in these cases the direction of shift ik also invariable regardless of the type 
of earth model. Israel et al. (1973) used a different method to obtain the dis- 
placement of the pole for model M 3 , in general the results of this investigation 
are in agreement with those given by them. 
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Table 11.1 


Pole Shift (pm) 


20 km 60 km 100 km 

200 km 


Vertical 
Strike- Slip 

' V 

7.3 

17.4 

18.3 

18.4 

m 3 

8.1 

17. ,3 

16.. 6 

17.8 

Vertical 

P 

2.0 


13.1 

27.7 

Dip-Slip 

U 3 

2.1 

7.7 

13.8 

28.9 

Dip -Slip 

P 

19.4 

45.3 

46.2 

45.4 

on 

45° Plane 

•m 3 

22 , 6 

4g. 7 

43.. 7 

44,3 


Table 11.2 


Angle (Degrees) 


20 km 60 km 

100 km 20 pkm 


Vertical 

P 

-147, 6 

-147, § 

-147, 6 

-14-7. j 6 

Strike-Slip 

m 3 

-147.6 

-147 . § 

-147.6 

> v- 

-147.6 

Vertical 

p 

91,0 

91*0 

91,0 ~ 

91. P 

Dip-Slip 


91.0 

91,0 

91.0 

** - Y ' 

pi, P 

Dip-Slip 

on 

p 

-35.6 

-38.0 

-43.2 

'“52,2 

4.5° Plane 

M 3 ' 

-28.4 

.-33,0 

-35.7 

^ i P 

-47,1 
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Table 11.3 
ALOD (see) x 10 8 




20 km 

60 km 

100 km 

200 km 

Vertical 

P 

127 

302 

317 

319 

Strike-Slip 

m 3 

140 

300 

287 

308 

Vertical 

P 

-52 

-189 

-339 

-718 

Dip-Slip 

m 3 

-56 

-199 

-357 

-748 

Dip-Slip 

on 

' 45° Plane 

p 

741 

1686 

1627 

1404 

M 3 

910 

1815 

1662 

1482 



Table 11.4 





> 

O 

to 

O 

xlO 10 





20 km 

60 km 

100 km 

200 km 

Vertical 

p 

-0.073 

-0.174 

-0.182 

-0.184 

Strike-Slip 

m 3 

-0.080 

-0.172 

-0.165 

-0.177 

Vertical 

p 

0,030 

0.109 

0.195 

0.413 

Dip-Slip 

m 3 

0.032 

0.114 

0.205 

0.430 

Dip-Slip 

on 

45° Plane 

p 

-0.426 

-0.971 

. -0.937 

-0.808 

M 3 

-0.524 

-1.045 

-0.957 

-0.853 
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Table 11, 5 
AC 2 j x 10 10 


20 km 

60 km 

100 km 

,20.0 km 


Vertical 

P 

-0,075 

.-.0.177 

-0.186 

-0.187 

I 

Strike-Slip 

m 3 

-0.082 

-0.17.6 

-0.168 

-0.180 

Vertical 

P 

-0.000£ 

-0,001 

-0.002 

-0. 00.6 

Dip-Slip - 

M 3 

-0.0004 

^0.001 

-0.003 

-0.0.06 

Dip- Slip 

rin 

P 

.0.190 

0.429 

0. 105 

0.334 

Ull 

45° Plane 

m 3 

0.239 

0.4.70 

0.427 

0.362 


Table 11. .6 


as 21 x lO' i0 


20;km ,60 km 100-km 

20O ! km 


Vertical 

P 

^0.047 

-=0.112 

-OL M8 

-0.119 

Strike-Slip 

m 3 

-0. 052 

— 0. I'll 

-0.:107 

-0.114 

Vertical 

P 

ft. 024 

Q. 088 

0.157 

0.333 

Dip-Slip 

M 3 

0.026 

0.092 

0.166 

,0.347 

Dip-Slip 

AVI ^ 

p 

-0.>13,6 

-0.336 

-0.380 

-0.431 

on 

45° Plane 


-0.4.29 

-0,306 

-0.307 

-0.390 
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Table 11.7 
•AC 22 x 10 11 




20 km 

60 km 

100 km 

200 km 

Vertical 

P 

-0.394 

-0. 933 

-0.979 

-0.987 

Strike-Slip 

M 3 

-0.432 

-0. 928 

-0.888 

-0.952 

Vertical 

P 

0.030 

0.110 

0.198 

0.418 

Dip-Slip 

m 3 

0.033 

0.116 

0,208 

0.436 

Dip-Slip 

on 

45° Plane 

p 

0.283 

0.685 

0.750 

0.808 

m 3 

0.289 

0.650 

0.639 

0.751 



Table 11. 8 
AS 22 x 10" 





20km 

60 km 

100 km 

200 km 

Vertical 

p 

-0,839 

-1. 984 

-2.081 

-2.099 

Strike-Slip- 

M 3 

-O'. 919 

-1.972 

-1.887 

-2.024 

Vertical 

p 

-0.070 

-0.253 

-0.453 

. -0.958 

Dip-Slip 

m 3 

-0.075 

-0.265 

-0.477 

-0.998 

Dip-Slip 

p 

0.603 

1.457 

1.595 

1.719 

on 

45° Plane 

M 3 

0.614 

1.383 

1.358 

1.597 
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FIGURE CAPTIONS 
Figure 10.1. y, vs. r. Model P 

Figure 10.2. y 2 vs. r. Model P 

Figure 10.3. y 3 vs. r. Model P 

Figure 10.4. y 4 vs. r, Model P 
Figure 10. 5. y 5 vs. r. Model P 
Figure 10, 6. y 6 vs. r, Model P 
Figure 10. 7. vs. r, Model M 3 
Figure 10.8. y 2 vs. r. Model M 3 - 
Figure 10. 9. y 3 vs. r. Model M 3 
Figure 10,10. y 4 vs. r, Model M 3 
Figure 10. 11. y 5 vs. r, Model M 3 


Figure 10.12. y 6 vs. r, Model M 3 



